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ABSTRACT 

Several mechanisms have been proposed to explain the alignment of grains 
with the interstellar magnetic field, including paramagnetic dissipation, radiative 
torques, and supersonic gas-grain streaming. These must compete with disalign- 
ing processes, including randomly directed torques arising from collisions with 
gas atoms. I describe a novel disalignment mechanism for grains that have a 
time-varying electric dipole moment and that drift across the magnetic field. 
Depending on the drift speed, this mechanism may yield a much shorter dis- 
alignment timescale than that associated with random gas atom impacts. For 
suprathermally rotating grains, the new disaligning process may be more po- 
tent for carbonaceous dust than for silicate dust. This could result in efficient 
alignment for silicate grains but poor alignment for carbonaceous grains. 

Subject headings: ISM: dust 

1. Introduction 

The phenomenon of starlight polarization was discovered over fifty years ago (Hall 1949; 
Hall & Mikesell 1949; Hiltner 1949a, 1949b) and was quickly attributed to dichroic extinction 
by a population of non-spherical, aligned interstellar dust grains. In the ensuing decades, 
many detailed investigations of potential alignment mechanisms were conducted, but to date 
no successful theory has been fully elaborated. See Lazarian (2003) and Roberge (2004) for 
recent reviews of grain alignment. 

Davis & Greenstein (1951) developed a promising model based on dissipation resulting 
from the rotation of a paramagnetic grain through the static interstellar magnetic field B. 
Since Fe is abundant and heavily depleted in the interstellar medium (ISM), it is likely that 
at least some grains are paramagnetic. Davis & Greenstein assumed that the grain rotation 
is excited by elastic impacts with gas atoms. In this case, the energy in rotation about any 
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grain axis is ~ \kT gsS) where k is Boltzmann's constant and T gas is the gas temperature; 
thus, this motion is called "thermal rotation" . 

In addition to exciting grain rotation, the random torques imparted by gas atom impacts 
disalign grains. Suppose each impact delivers angular momentum mva, where m is the mass 
of the gas atom, v is its typical speed, and a is the grain size (which we will take to be the 
radius of a sphere of equal volume). After iV collisions, the angular momentum axis will 
have suffered an angular displacement A9 ~ [mva/ J)N 1 / 2 ; J is the grain's total angular 
momentum. For thermal rotation, mv 2 /2 ~ J 2 /2I, where / ~ Ma 2 is the grain's moment 
of inertia (M is the grain's mass), implying A9 ~ (iVm/M) 1 / 2 . For A9 ~ 1 rad, N ~ M/m. 
Thus, the timescale for disalignment, Tdi SiCO i, equals the time for the grain to collide with its 
own mass in gas. If the gas has number density n and the grain, with mass density p, does 
not drift relative to the gas, then 

The disalignment timescale is expected to be shorter than the Davis-Greenstein alignment 
timescale, unless the grains have superparamagnetic inclusions (Jones & Spitzer 1967). 

Thus, other alignment mechanisms were sought. Martin (1971) pointed out that a 
charged, spinning grain has a magnetic dipole moment /x. Dolginov & Mytrophanov (1976) 
showed that the Barnett effect (i.e., the tendency for a spinning paramagnetic solid to acquire 
a magnetization parallel or anti-parallel to its angular velocity to) can provide a much larger 
moment. Specifically, the Barnett magnetic moment /^Bar = Xo^V/lg, where xo is the static 
magnetic susceptibility, 7 9 is the gyromagnetic ratio of the microscopic magnetic dipoles 
that are responsible for the grain's paramagnetism, and V is the grain volume. The torque 
= fix B causes the grain to precess about the interstellar magnetic field B. The precession 
rate is fast enough that, even if the actual grain alignment mechanism did not involve B, 
the observed polarization would be either parallel or perpendicular to B. 

Purcell (1979) noted that the Barnett effect also underlies a mechanism for dissipating 
grain rotational energy. For fixed angular momentum J, there are many rotational states 
available to a rigid body, with different rotational energies. If the grain is not executing 
steady rotation about one of its principal axes, then ui varies periodically in grain body 
coordinates. The magnetization lags u> and the resulting paramagnetic dissipation, which 
Purcell called "Barnett dissipation", transfers rotational energy to the thermal reservoir 
provided by the grain vibrational modes. This drives the principal axis of greatest moment 
of inertia, d\, into alignment with J, since this configuration has minimum rotational energy 
for given J. Note that this mechanism acts even in the absence of an external magnetic field. 
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Purcell (1975, 1979) also noted that grains are subject to systematic torques fixed 
in grain body coordinates and that these torques can spin the grains up to suprathermal 
rotational speeds. The most important torque identified by Purcell results from the formation 
of H 2 at special surface sites, with subsequent ejection from the grain. This discovery restored 
the appeal of the Davis- Greenstein mechanism, since suprathermally rotating grains are 
largely impervious to disalignment by random gas atom impacts. However, the distribution 
of H 2 formation sites on the grain surface is expected to change (e.g., due to accretion of 
material from the gas) on a timescale short compared with the alignment time. Consequently, 
the associated torque will sometimes be oriented antiparallel to J, spinning the grain down 
to the thermal rotation regime. Spitzer & McGlynn (1979) found that grains are effectively 
disaligned during these "crossover" episodes. 

Lazarian & Draine (1997; 1999a, b) reexamined the physics of crossovers, including the 
effects of grain thermal fluctuations. Thermal fluctuations, in which energy is spontaneously 
transferred from the vibrational modes in the solid to grain rotation (at constant J), are 
associated with the internal dissipation mechanisms (which we assume to be dominated by 
Barnett dissipation). 

Before proceeding, we must first summarize some key features of the torque-free rota- 
tional motion of asymmetric rigid bodies; see §2.5 of Weingartner & Draine (2003, here- 
after WD03) for a more detailed description. Suppose the inertia tensor has eigenvalues 
I\ > I2 > -^3) with principal axes d±, a 2 , and a 3 . It is convenient to express the rotational 
energy in terms of the following dimensionless parameter: 

q =— ■ (2) 

Note that 1 < q < I1/I3. The grain rotational state is not completely specified by (J, q). For 
any (J, q), the grain is in one of two possible "flip states". These flip states are defined with 
respect to one of the principal axes; which one depends on the value of q. The grain is in 
the positive (negative) flip state with respect to &i if the component of the angular velocity 
along that axis, oui, remains positive (negative) throughout the grain rotational motion. 

Lazarian & Draine (1999a) identified an important consequence of thermal fluctuations, 
which they called "thermal flipping". Suppose a grain initially has q < I1/I2 and is in 
the positive flip state with respect to di. A thermal fluctuation may add enough energy 
to the rotational motion so that q > I1/I2 (in which case the grain will be in one of the 
flip states with respect to 03). If the evolution subsequently returns q to a value < h/h, 
then the grain may end up in either the positive or negative flip state with respect to d±] 
in the latter case a thermal flip has occurred. When a grain flips, any vector fixed in 
grain-body coordinates, like the H 2 formation torque, reverses direction in space. Thermal 
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fluctuations become more pronounced as the grain angular momentum J decreases. Lazarian 
& Draine (1999a, b) found that thermally rotating grains with a < 1/im undergo such rapid 
flipping that the H 2 formation torque may never have the opportunity to spin a grain up 
to suprathermal rotation, a situation that they called "thermal trapping". Thus, it appears 
that Davis-Greenstein alignment with H2 formation torque-induced suprathermal rotation is 
not viable. 

Draine & Weingartner (1996) found that radiative torques (due to the absorption and 
scattering of starlight by an irregularly shaped grain) can also yield suprathermal rotation. If 
the interstellar radiation field is modestly anisotropic, then (1) the radiative torque does not 
reverse direction on each grain flip, perhaps overcoming the problem of thermal trapping, and 
(2) radiative torques can directly align grains, on a shorter timescale than that of the Davis- 
Greenstein mechanism (Draine & Weingartner 1997). Thus, the radiative torque alignment 
scenario can be summarized as follows. The radiative torque, averaged over grain rotation 
and precession, drives the grain to suprathermal rotation and aligns J with respect to B. 
On a much faster timescale, Barnett dissipation aligns J with d\. 

Suppose a grain has an electric dipole moment p, in addition to its magnetic dipole 
moment, and that it drifts with velocity v relative to the gas and magnetic field. This grain 
will experience a torque 

r p = p x Q x b) , (3) 

where c is the speed of light. Consider a coordinate system x, y, z in which B = Bz and 
v = v\\z + v±y. The field axis z also serves as the polar axis for spherical coordinates, with 
polar angle 9 and with the azimuthal angle <f> measured relative to x. Suppose in addition 
that di, n, and p lie parallel or anti-parallel to J: fi = jijJ and p = pjj; both jij and 
pj can be positive or negative. Then, adding the torques due to the magnetic and electric 
dipoles yields 

T = fijBJ x [z + Tx] (4) 

with 

T = ^ . (5) 

HjC 

Thus, the grain precesses not about z, but about another axis tilted at polar angle 5 = 
tan -1 |T|. The azimuthal angle of this axis is <fi = (ir) when T > (T < 0). Also, the 
precession rate is increased by a factor (1 + Y 2 ) 1 / 2 . 

In general, and p are not parallel or anti-parallel to J. How does this affect the 
dynamics? The torques should be averaged over grain rotation, since this occurs on a much 
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faster timescale than any other relevant process. Denoting such averages with a bar, 

f -. = **B=(|>xB=(£)£xB (6) 

(WD03, §3.2). If the Barnett effect is responsible for the magnetic dipole moment, then 
(fJ-fio) = XoV/lg- Note that we have ignored the component of the magnetization that is out 
of phase with u>; this is justified in §3.2 of WD03. Note also that f M is a factor q larger than 
its value when J || d\ and it does not depend on the flip state. 



Similarly, 

From §3.4 of WD03, 
where 



Mq) = 



f p = p x (X x B) . (7) 
p(q,±) =±(p-a i )f i (q)J (8) 

/*-/3?n 1/a 



h-h 



(9) 



2F(tt/2, kf) 



the ± refers to the flip state, % — 1 (3) when q < h/h (q > h/h), F is the elliptic integral 
of the first kind (see eq. 15 of WD03), and 



(J 2 -/ 3 )(g-l) 



K ~ (h-m-hq/h) (10) 

when q < I1/I2 and its inverse when q > 

Thus, if the grain rotational state remains fixed, then the dynamics is identical to that 
for the case that J, /x, and p all lie along d\ (eqs. 4 and 5), with 

and 

p J = ±(p-d t )f l (q) . (12) 
Figure 1 displays fi(q)/q for a grain with h : I 2 : h — 1 : 0.8 : 0.7. 

The role of the electric dipole is insignificant (i.e., T <^ 1) when the grain rotates 
suprathermally and the drift speed perpendicular to the magnetic field is low. Both of 
these conditions have traditionally been assumed, since grain alignment seemed to require 
suprathermal rotation and gyrorotation and drag were thought to suppress drift across field 
lines. 



However, two recent studies call these assumptions into question. First, WD03 gener- 
alized the Draine & Weingartner (1997) analysis of grain alignment by radiative torques. 
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Motivated by the fact that radiative torques drive suprathermal rotation, the earlier treat- 
ment adopted the assumption that J || a±. However, radiative torques also naturally drive 
grains through crossovers; the importance of thermal fluctuations and flips during these 
episodes invalidates this assumption. Thus, WD03 dropped this requirement. Surprisingly, 
they found that radiative torques can align J with respect to B without driving the grain 
to suprathermal rotation. For purposes of computational speed, WD03 did not adopt a 
realistic spectrum for the interstellar radiation field. As a result, they could not estimate 
the timescale for alignment in the thermal rotation regime. Thus, it is not yet clear whether 
the timescale for disalignment by random gas atom collisions is longer or shorter than the 
alignment time in this regime. 

Second, Yan, Lazarian, & Draine (2004) found that interstellar magnetohydrodynamic 
turbulence can accelerate grains to speeds as high as ~ lkms -1 across magnetic field lines 
in the cold neutral medium. Thus, mildly supersonic motion is possible. In this case, the 
angular momentum imparted by collisions with gas atoms tends to lie perpendicular to the 
drift velocity. Yan et al. (2004) suggested that the associated "mechanical alignment" (Gold 
1951, 1952) may be more important than previously thought. Mildly supersonic drift also 
yields mildly suprathermal rotation. 

2. The Role of the Electric Dipole in Disalignment 

Suppose pj (eq. 12) is constant. In this case, the only consequence of the electric dipole is 
to tilt the precession axis relative to the magnetic field direction. Because of gyrorotation, the 
precession axes of individual grains would likely be uniformly distributed in azimuth (about 
the field direction). Thus, the observed polarization is still either parallel or perpendicular 
to the field, but is diluted. Consider the simple case where the magnetic field lies in the 
plane of the sky and the grains are aligned such that the precession angle (with respect to 
its instantaneous precession axis) is zero. For T C 1, the polarization is hardly affected. 
As T increases, the dilution increases as well. When T = 1, the precession axis is tilted 
45° from the field direction. Two grains that are 180° out of phase in their gyrorotation 
produce orthogonal polarization vectors. Thus, the polarization is completely suppressed 
when T = 1. As T increases beyond 1, the dilution decreases, but does return to zero as 
T — > oo. This polarization dilution is identical to that which occurs when T = and the 
grains are aligned with a non-zero precession angle (with respect to the field direction). 

In the above paragraph, we assumed that pj is constant; in fact, it exhibits variations 
on a relatively short timescale. Whenever the grain experiences a discrete charging event 
(e.g., the capture of an electron from the gas or the photoejection of an electron), p will 
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change. It is not clear how the orientation of p in grain body coordinates is determined. It 
may depend primarily on the overall grain geometry. Alternatively, it may be sensitive to 
localized irregularities where charge can be concentrated. In the former case, the orientation 
of p is probably fairly constant. In the latter case, it may change with each discrete charging 
event, and may even reverse direction with respect to d\ from time to time. Finally, p may 
harldy vary at all (in grain body coordinates). This would be the case if it originates primarily 
in the random distribution of polar constituents rather than the asymmetric distribution of 
excess charge. 

Even if p is fixed in grain body coordinates, p,j will still vary, if the rotational state (i.e., 
q and/or the flip state) varies (eq. 12). The value of q can be affected by both the action of 
external torques and internal relaxation processes, like Barnett dissipation and fluctuations. 
Of course, thermal flipping yields changes in the flip state. 

The dynamical equation for q is 



(WD03, §2.4); internal relaxation is manifested in the J • doj/dt term. The electric dipole 
does not affect q, since F p ■ u) = and J ■ T p = 0. This would also be true for the magnetic 
torque if we ignore the lagging component of magnetization, as in equation (6). The lag- 
ging component does yield a non-vanishing contribution, but it is generally small compared 
with the Barnett dissipation rate for the grain sizes considered here (§3.2 of WD03). The 
contribution of other torques, like the radiative torque, are only significant on much longer 
timescales and can be ignored. Thus, variations in q are dominated by Barnett dissipation 
and fluctuations, characterized by a timescale rear- Of course, the thermal flipping timescale 
r t f (i.e., the inverse of the probability per unit time that a flip occurs) must be related to 

fBar- 

Variations in /ij and p j (and hence T) cause the precession axis to vary, resulting in 
a potent mechanism for disaligning grains. This process is most simply demonstrated if we 
assume that the magnitudes of \ij andpj are constant, but that the grain undergoes thermal 
flipping, reversing the sign of pj. Suppose that initially the grain is in the positive flip state 
with respect to d\ and J is precessing, with polar angle $i, about its axis, tilted at angle 5 
relative to the magnetic field direction. Following a thermal flip, the precession axis abruptly 
changes to the one associated with the negative flip state. The angle between these two axes 
is 26. The grain precesses about this new axis for some time and then undergoes a thermal 
flip. At that point, it resumes precessing about the original axis, but the polar angle now 
will not equal Q\. After many flips, the grain's precession angle may differ substantially from 
its initial value. 




(13) 
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In order to estimate the disalignment timescale, it is convenient to write the equations 
of motion for J in the coordinate system (9, 0) with the magnetic field as the polar axis: 

^ = Qo[l - Tcot#cos(0 + cj gyro t)] (14) 
d9 

— = -fi o Tsin(0 + ^ gyro t) , (15) 

where Qo = —f^jB/J is the precession angular velocity in the absence of an electric dipole 
and o;g yr o is the gyrofrequency. In the limit r t f — > 0, the terms proportional to T average to 
zero and the grain simply precesses about the magnetic field direction; 9 remains fixed. The 
timescale for gyrorotation is given by 

u - ~ 24 * 102 fey (o^v)" 1 (^r yr • < ie) 

where U is the grain potential. For most cases of interest, this is much longer than the 
timescales for precession and flipping. 

Suppose r t f is much shorter than the precession timescale, ily 1 = |fio| _1 (l + T 2 ) -1 / 2 . 
Consider the evolution over a time interval dt\ + rft 2 , during which the grain spends time dti 
in the positive flip state and dt2 in the negative flip state. From equation (15), the resulting 
deviation in polar angle d9 ~ |f2o| T \dti — dt 2 \ ~ |fio| Y r tf- After iV double flips, the angular 
deviation A# ~ N 1 ^ 2 d9. For A# ~ 1 rad, ~ d9~ 2 . Thus, the disalignment timescale 

T^tf ~ Nr ti ~ r-^r 2 ^ 1 , if r ti <^\n \-\i + T 2 )' 1/2 ■ (17) 

In the opposite limit, where r t f ^> Q^ 1 , the grain executes many full precession cycles, plus 
one partial cycle, before flipping. Since the deviation in 9 over a full precession cycle equals 
zero, d9 is due to the difference in the final, partial precession cycles for the two flip states. 
Thus, d9 ~ |Q | T fix 1 = T(l + T 2 )- 1 / 2 ; this implies 

r d is,tf ~(1 + T^ 2 )r tf , if r tf »|fi r 1 (l + T 2 )- 1 / 2 . (18) 

Note the opposite dependence of the disalignment timescale on the flipping timescale in the 
above two regimes; Tdi s ,tf oc r t f (t^ 1 ) when flipping is slow (fast) compared with precession. 

Of course, this scenario is simplified. There are not sudden transitions between two 
values of T, with the same magnitude but opposite signs. Rather, the grain experiences 
smaller variations in T on a shorter timescale, rear- Flipping is a cumulative result of 
these variations. We do not expect this complication to affect the above order-of-magnitude 
estimates for r^tf- As long as the flip state remains the same, the sign of d6 will also remain 
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the same (in the limit that r t f <C Qy 1 ), even as T varies. Thus, it is still appropriate to 
consider the evolution between two flips to constitute a single step in the random walk. The 
only modification to the disalignment timescale is that T should be replaced with (T), its 
thermal average over a distribution of rotational states (i.e., values of q) in a single flip state 
about a\. For a grain with temperature Td, Ii : I 2 : h = 1 : 0.8 : 0.7, and J 2 / \2IikTd) = 6, 
(T)/T(g = 1) f« 0.7; T reaches its maximum value when q — 1. Thus, this modification is 
not very important. (See §4.1 of WD03 for a derivation of the density of states as a function 
of q, which is needed to perform the thermal average.) 

Although a detailed model of Barnett relaxation and thermal flipping has not yet been 
developed, Lazarian & Draine (1999a) suggested that 
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where 



In the case of thermal rotation, (J/J ) 2 ~ T gas /T d ; thus r tf is longer than r Bar by one to two 
orders of magnitude for thermally rotating grains in the cold neutral medium. 

If discrete charging events substantially alter the orientation of p in grain body coordi- 
nates, then this process may also yield a form of flipping. In this case, p flips in grain body 
coordinates rather than the grain orientation flipping with respect to J. The consequence 
for grain disorientation is the same. If Td c is the timescale on which discrete charging events 
occur and the fractional change in pj on each such event is ~ Z^ 1 (Z is the excess number 
of elementary charges on the grain), then the associated "flipping" timescale Tf )dc ~ Z 2 r dc . 
If r f,dc < r ti, then Tf 5 dc should replace r t f in the estimate of the disalignment time. 

Thus, the disalignment timescale depends on (1) the drift speed v± across the magnetic 
field, (2) the precession timescale \tto\, (3) /ij and pj, and (4) r t f (which is related to rear) 
or rf ; dc- In the next section, we will estimate the values of these quantities for interstellar 
grains. 



3. Dipole Moments and Relaxation Timescales 

The magnetic dipole moment of a uniformly charged sphere with radius R, total charge 
q, and rotational speed u is /i c harge = R 2 quj/3c. Thus, 

Charge w 1.1 x 10~ 24 ( ( j^^J (l^l) ^ m statcoulomb cm , (21) 
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where U is the grain potential and K m accounts for deviations from spherical shape. Wein- 
gartner & Draine (2001) estimated that U ~ 0.3 V for carbonaceous and silicate grains with 
a ~ 0.1/im in the cold neutral medium (see their fig. 10). The thermal rotation rate for a 
sphere with radius a is given by 

V 8vrpa 5 J V 3 S cm / V 100K / V - 1 ^ 111 / 

where p is the density of the grain material. 

The Barnett magnetic moment /iBar = Xo^V/lg- The susceptibility is given by Curie's 
Law (see, e.g., Morrish 1980): 



n 7 2 ft 2 J( J + 1) 
SkT d 



Xo = ^77F " > (23) 



where J is the angular momentum quantum number of the paramagnetic ion or nucleus, n 
is its number density, and is the dust temperature. 

First, consider silicate grains with structural unit MgFeSi04. The magnetization is 
dominated by the Fe ions. Thus, we take J = 5/2 and 7 9 = —1.76 x 10 7 s _1 G _1 , yielding 
Xo ~ 5 x 10 _3 (T rf /15 K)~ x (Draine 1996). For silicates with less Fe incorporation, xo would 
be lower. With these estimates, 

|" B „|( S il) " L2 X 1<r " (wk) 1 (ol^)' (wpd rtatcoulombcm . (24) 

Next, consider carbonaceous grains with H incorporation at the ~ 10% level (i.e., the 
number density of H is « 10 22 cm~ 3 ). The paramagnetism is dominated by the H nuclei, 
with 7 9 = 2.675 x lO^G" 1 and J = 1/2, implying X o ~ 9.6 x 10- 10 (T d / 15 K) _1 . Thus, 
we estimate 

/Wcarb) « 1.5 x KT 23 (^-^ statcoulomb cm . (25) 

If the carbonaceous grains are not hydrogenated to any significant extent, then the nuclear 
paramagnetism is contributed instead by 13 C, which is present at the ~ 1% level and has 
7 9 = 6.73 x 10 3 s _1 G _1 and J = 1/2. In this case, /i Bar would be reduced by a factor 
n( 1 H)7 s ( 1 H)/n( 13 C)7 9 ( 13 C) = 40 compared with the above estimate. On the other hand, 
dangling bonds in hydrogenated amorphous carbon may yield electron paramagnetism with 
a spin density ~ 10 19 g _1 (Esquinazi & Holme 2005). This is a factor 10 3 less than the Fe 
density in the model silicate considered above. Thus, the Barnett moment associated with 
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these dangling bonds may exceed the estimate for hydrogenated carbon in equation (25) by 
an order of magnitude, depending on the value of j g . 

If the electric dipole moment is due to an asymmetric distribution of excess charge on 
the grain (characterized by an uknown parameter K e ), then 

p = 9 a Kc « 1.0 x 10- (J^) (gJL-)' (J^) statcoulombcm . (26) 

We have adopted a fiducial value of n e 10~ 2 . This seems conservative, since the number 
of excess charges on a grain with U = 0.3 V and a = 0.1/zm is only 20. In addition, if the 
grain is composed of N randomly arranged polar constituents, each with dipole moment p , 
then we would expect an additional contribution to p ~ A rl / 2 p - This may, in fact, dominate 
the contribution due to excess charge (see, e.g., eq. 11 in Draine & Lazarian 1998). 

Observations of the wavelength dependence of starlight polarization reveal that rela- 
tively large grains (a > 0.1/xm) are aligned while smaller grains (a < 0.05/xm) are not (e.g., 
Kim & Martin 1995). For aligned grains, the Barnett magnetic moment dominates that due 
to the spinning charge (compare eqs. 21 and 25). For hydrogenated carbonaceous grains, the 
above estimates yield 



/iBarC yO.lfimJ VlO^- 1 / \ 15K / V0-3Vy V10- 2 / Vlkms" 

Assuming thermal rotation (with T gas = 100 K) and a = 0.2/zm (and that T, U, K e , and 
v± have the values in eq. 27), pv±/ /iBarC ~ 370 for hydrogenated carbonaceous grains and 
~ 0.046 for silicate grains. Of course, these estimates are highly uncertain. 

The precession rate 

~ ~T~ ~ °'° 31 (lO- 23 statcoulombcm) \ 57Ig) V3(?cm" 3 ) (o.lum) (lO^- 1 ) 



yr 
(28) 



Since \i oc u, f2 is independent of u>. 

We adopt the order-of-magnitude estimate of the Barnett relaxation timescale from 

V 



WD03 (their eq. 42). With the additional assumptions that Iy w 2I 3 and J 3 ~ tpVa 2 , 



i 5T 2 xo 

where 7i and T 2 are the spin-lattice and spin-spin relaxation times, respectively. With 
equation (23), this yields 

TbL ~ 7 x 10 3 ( 1Q22 cm _ 3 ) (lo - ^) (ol^i) (i5k) ( 



3 gem" 



)■' 


"J(J+ 1)" 
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where co>5 = c<j/10 5 s -1 . Note that rear oc T 2 but does not depend on j g . Thus, the constituent 
that dominates the Barnett relaxation may be different than that which dominates //Bar- The 
fiducial values in equation (30) are all intended for hydrogenated carbonaceous grains, for 
which Lazarian & Draine (1999b) estimated 7\ ~ T 2 ~ 10~ 4 s. For non-hydrogenated car- 
bonaceous grains, Tear will be longer. If silicate grains are hydrogenated, then equation (30) 
is applicable to them as well. Otherwise, the 29 Si nucleus (with 7 9 = —5.32 x 10 3 s _1 G _1 , 
J = 1/2, and n ^ 5 x 10 20 cm" 3 ) will likely dominate Barnett relaxation; this would increase 
r Bar by a factor ps 20. For thermally rotating grains, the contribution of nuclear paramag- 
netism to Barnett relaxation dominates that of electron paramagnetism, since T 2 is orders 
of magnitude shorter for electron paramagnetism. 

The discrete charging timescale can be simply approximated as the inverse of the rate 
at which electrons collide with the grain: 

(31) 

where x e is the electron fraction, n e /nn, and the final term accounts for Coulomb focusing. 
We also estimate the timescale for the flipping of p in grain body coordinates (if discrete 
charging events can yield substantial changes in p) as r^dc ~ Z 2 t& c with Z ~ 20(a/0.1/xm) 
for a ~ 0.1/zm. 



4. Consequences 



Consider hydrogenated carbonaceous and silicate dust grains with U ~ 0.3 V in the cold 
neutral medium. For the carbonaceous dust, we estimate (T) ~ 44(a/0.1/zm) L5 (c<j/c<> r )~ 1 
(i.e., 1/3 of the value estimated for pv±/fiB aT c in §3). From equations (19), (30), (28), and 
(31), and assuming the fiducial values therein, 
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(32) 



(33) 



-l 



21(a/0.1/zm) 2 yr and r fidc ~ 6 x 10 3 yr. Similarly, for silicate dust (T) 



\Q 

0.0056(a/0.1/im) 1 - 5 (cj/^ T ) _1 and |^ | _1 ~ 0.0027(a/0.1/W) 2 yr. The contribution of elec- 
tron paramagnetism to can be ignored since it is negligible for low to moderate to /cut- 
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For high u/ljt, where electron paramagnetism may be important for Barnett dissipation, 
the thermal flipping timescale is too long for the associated disalignment to be significant. 
Thus, we adopt equations (32) and (33) for silicate grains as well as for carbonaceous grains. 

We employ equation (17) to estimate the timescale for disalignment associated with 
thermal flipping, Tdi S ,tf, when r t f < Ifiol^l+Y 2 ) -1 / 2 and equation (18) when r t f > |Qo| _1 (l+ 
Y 2 ) -1 / 2 ; similarly for the discrete charging disalignment timescale, Tdi S ,dc- Figures 2 and 3 
show rd is t f and Tdi S ,dc versus lv/out, for grains with a = 0.1/xm and 0.2/im, respectively. For 
the lowest values of u, Tdi Si tf is identical for carbonaceous and silicate grains. In this regime, 
r t f < fig *(1 + Y 2 )^ 1 / 2 for both grain types. Although Y and f2o are very different for the 
two compositions, the disalignment timescale depends on the product Yf2 , in which the 
magnetic dipole moment cancels. 

For grains with a = 0.2/im, the precession timescale always exceeds Tf ^c, for both car- 
bonaceous and silicate grains. Thus, the discrete charging disalignment timescale is identical 
for the two grain types. In contrast, Tdi S ,dc is a factor 5 longer for silicate dust than for 
carbonaceous dust when a = 0.1/im. For silicate grains, r^dc exceeds the precession time; 
the opposite is true for carbonaceous grains. Thus, Tdi S ,dc = Y " 2 Tf i d c (^/^T) 2 for silicate 
grains and Tdi S ,dc = Y 2 Qq 2 t^\ c {uj/ujt) 2 for carbonaceous grains; Y = T(cu = ujt)- 

Suppose discrete charging events do not substantially change the orientation of p in grain 
body coordinates. If the adopted parameter values are accurate, then thermally rotating 
grains of both compositions are disaligned extremely rapidly. The scenario of WD03, in which 
radiative torques align grains with uj u>t, would seem highly unlikely. Alignment with 
suprathermal rotation may work, but disalignment associated with the electric dipole would 
play a major role in the crossover dynamics. Since Tdi s ,tf oc (pv±)~ 2 , errors in the assumed 
values of p and v± could invalidate these conclusions. Although the estimated value of p 
seems conservative, a detailed model of the electric dipole is needed. Also, confirmation that 
magnetohydrodynamic turbulence can accelerate grains to v± ~ lkms -1 (Yan et al. 2004) 
is needed. Finally, a detailed model of thermal flipping is needed. If the steep increase in 
r t f is shifted to somewhat lower values of u/uj t , then alignment of thermally rotating grains 
may be tenable, even if pv± is large. 

If discrete charging events can lead to flips of p in grain body coordinates, then even 
suprathermally rotating grains may be subject to significant disalignment. The Davis- 
Greenstein alignment time > 10 6 yr for a > 0.1/im (see, e.g., eq. 71 in WD03). Draine & 
Weingartner (1997) found that radiative torques yield aligned grain states with lo/uo t ~ 100 
and alignment timescales > 10 5 yr. If the adopted estimate for r^dc is accurate, then Davis- 
Greenstein alignment with uj/ujt < 100 seems at best marginal, even if crossovers could be 
avoided. Alignment by radiative torques fares better, but even in this case disalignment may 
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play a significant role. 

Suppose the adopted value of the flipping timescale due to discrete charging, r^dc, is 
too short. When a = 0.1/xm, higher values of r fjdc would yield higher values of the disalign- 
ment timescale Tdi S ,dc for silicate grains (eq. 17) and lower values of r^dc for carbonaceous 
grains (eq. 18). If r fidc were sufficiently large that it exceeds the precession timescale for 
silicate grains with a = 0.2/im (but not for carbonaceous grains), then disalignment would 
occur more rapidly for carbonaceous grains than for silicate grains of this size. If Tf ;( jc were 
larger than the adopted value by a factor of ~ 40, then 7di Si dc <^ 10 7 yr for silicate grains, 
but Tdis,dc ^ 10 4 yr for carbonaceous grains, when a m 0.2yum and uj/uj^ m 10 2 . Thus, 
alignment may proceed more quickly than disalignment for silicate grains, but vice versa 
for carbonaceous grains. This could explain the observation that silicate grains are aligned 
while carbonaceous grains may not be (Whittet 2004). Again, a detailed model of the electric 
dipole moment in interstellar grains is needed to clarify this possibility. 

With the above estimates, (T) ~ 1 for carbonaceous grains with a 0.2 fim and u/ut ~ 
100. Thus, even if these grains do undergo efficient alignment, the polarization may be 
significantly diluted. 

5. Summary 

1. When a grain with an electric dipole moment drifts across the magnetic field B, the grain 
precesses about an axis different from B. 

2. The importance of the electric dipole torque scales with the parameter T = pj v±/ fij c; 
T is expected to be larger for carbonaceous grains than for silicate grains. 

3. Because of gyrorotation, polarization still correlates with B when T > 0, but the polar- 
ization is diluted, with maximal dilution when T 1. 

4. Variation of the sign of T, due to thermal flipping or discrete charging, produces grain 
disalignment. The timescale for this disalignment process may be much shorter than that 
for disalignment by random collisions with gas atoms. 

5. In the case of suprathermal rotation, the disalignment timescale may be shorter for 
carbonaceous grains than for silicate grains. 

6. A detailed model of the electric dipole moment in interstellar grains is needed before firm 
estimates of the disalignment timescale can be made. 
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Fig. 1. — fi(q)/q, (see eq. 9), with i — 1 when g < -Z1//2 and i — 3 when g > I1/I2, for a 
grain with A : J 2 : J 3 = 1 : 0.8 : 0.7. 
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0.1 1 10 1 10 2 10 3 

Fig. 2. — Disalignment timescales for changes in p associated with thermal flipping (r^tf) 
and discrete charging (rdi S ,dc), for carbonaceous and silicate grains with a = 0.1/iin. The 
suprathermality is indicated by the ratio of the angular speed u to its thermal value, u>t- 



-19- 




0.1 1 10 1 10 2 10 3 



Fig. 3. — Same as fig 2, but for grains with a = 0.2//m. 



